Quantifying the Reversible Association of Thermosensitive Nanoparticles 
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Under many conditions, biomolecules and nanoparticles associate by means of attractive bonds, 
due to hydrophobic attraction. Extracting the microscopic association/dissociation rates from ex- 
perimental data is complicated by the dissociation events and by the sensitivity of the binding force 
to temperature (T). Here we introduce a theoretical model that combined with light scattering 
experiments allows us to quantify these rates and the reversible binding energy as a function of 
T. We apply this method to the reversible aggregation of thermoresponsive polystyrene/poly(N- 
isopropylacrylamide) core-shell nanoparticles, as a model system for biomolecules. We find that 
the binding energy changes sharply with T, and relate this remarkable switchable behavior to the 
hydrophobic-hydrophilic transition of the thermosensitive nanoparticles. 

PACS numbers: 



The association and self-organization of 
biomolecules pQ and nanoparticles [2) play a funda- 
mental role in many physiological processes [3], such 
as protein amyloid formation which is responsible for 
neuro-degenerative diseases [4], as well as in technologi- 
cal processes. In particular the assembly of nanoparticles 
has become a key step in the synthesis of nanomaterials 
with new optical and mechanical properties [5 ]. In 
order to understand and control all these association 
processes, it is essential to understand and control the 
association kinetics [H 0]. This is challenging in the 
case of biomolecules and nanoparticles where often the 
microscopic binding energy, which directly determines 
the rate of dimer formation, is comparable to the thermal 
energy ksT. This poses two major difficulties: (i) the 
experimental detection of the binding energy has to be 
accurate down to the fc^T scale, and (ii) the association 
kinetics is affected by dissociation events [5J [7J since the 
binding energy is comparable with the average kinetic 
energy of the particles/molecules. It has been recently 
shown that a clear understanding of protein aggregation 
under physiological conditions cannot be achieved 
without the understanding of aggregate dissociation [7J . 
A further complication arises from the fact that the 
binding energy can be very sensitive to changes of the 
solution parameters, such as T and pH [3] ■ 

In this Letter we propose a detection strategy that can 
overcome many of the difficulties to date in extracting 
quantitative information about the intrinsic rates of as- 
sociation and dissociation processes. The essence of the 
strategy is to combine a kinetic model with the analysis of 
dynamic light scattering (DLS) data of aggregation kinet- 
ics. Dynamic light scattering has proven a reliable tool 
for analyzing irreversible aggregation [T2], Here we show 
that it can be used for measuring reversible processes as 
well, but the data analysis requires a new model fully 
accounting for dissociation processes. Here we introduce 



such a model and demonstrate its applicability in light 
scattering experiments of thermosensitive nanoparticles 
in aqueous suspension. We also show how to quantify 
the T-dependent interaction energy between the parti- 
cles once the rates are extracted. 

Thermosensitive nanoparticles are interesting on their 
own and are the object of intense study because their 
size and interaction can be modified by changes of T [5] . 
They consist of a 52 nm radius solid polystyrene core 
onto which a polymeric network of crosslinked poly(N- 
isopropyl acrylamide) (PNIPAM) of T-dependent thick- 
ness (~ 50nm at T < 32°C and ~ 33nm at T > 32°C) 
is affixed (see Figjljl. At room temperature this ther- 
mosensitive shell is swollen by the dispersing agent wa- 
ter, which is expelled if the suspension is heated above 
the critical temperature T C =32°C of a volume transition. 
It is a well-established fact that these particles become 
attractive above T c [9| . This is due to the hydrophobic in- 
teraction between the PNIPAM-network for which water 
has become a poor solvent under these conditions. As a 
consequence of this, aggregation sets in which is entirely 
reversible upon cooling. 

The thermosensitive particles have been synthesized 
and characterized as described in [9]. Core-shell parti- 
cles with 5 mol.% crosslinking are used in this study. 
In salt-free solution these particles have a good stability 
at all temperatures. Adding salt, however, screens the 
electrostatic repulsion and induces aggregation above T c 
[9] . The reversible character of the aggregation is demon- 
strated in [8|. The kinetics was investigated by DLS us- 
ing a ALV 5000 light scattering goniometer (Peters) at 
a scattering angle of 90° and wavelength A = 632. 8nm. 
The experiment was done as follows: 2.3 mL of the latex 
solution (either at 2.72 x 10" 3 , 1.36 x 10~ 3 or 0.27 x 10~ 3 
wt.%) was equilibrated at the required temperature for 
20 min. Then 0.2 mL of a 0.625 M solution of KCl kept 
at the same T were quickly added to induce aggregation. 
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After homogenization the measurement was started and 
the average hydrodynamic radius Th of the colloidal con- 
stituents of the suspension was monitored as function of 
time (see Fig. 2). The time evolution of is linear in 
the early stage of association and for sufficient dilution 
and gives access to the association rate. If dissociation 
events occur on the time scale of observation, however, 
the measured association rate is an effective one (i.e. it 
contains the effect of dissociation and is therefore smaller 
than in the absence of dissociation) . In the following we 
develop a model which allows us to account for this effect 
and to extract the microscopic dissociation rate from this 
effective association rate from the DLS experiments. 
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FIG. 1: (color online). Structure of the thermosensitive par- 
ticles and of the effective square- well interaction, setting in at 
T > T c , as a function of the surface separation h. 

We start by considering the kinetics of reversible asso- 
ciation between two particles A to form a dimer A 2 



A + A ^ An 



(1) 



The association rate be denoted by k + and the dissocia- 
tion rate by fc_ . If we denote with n\ the concentration 
of monomers A at time t and with N the total concentra- 
tion of monomers at t = 0, the evolution of A is governed 
by the following equation 



dnAt) , , . 9 1 , 



\k-n x {t) (2) 



where we made use of the conservation condition: n 2 (t) = 
(N — m(t))/2, with n 2 the concentration of dimers A 2 . 
With the initial condition m(0) = N, Eq.(2) admits the 
following solution: 



2k, 



tanh(VU*/ 2 ) + B/VA 



1 + (B/VA) tanh(v0to/2) 



(3) 



2k + N. 



with A = k- (fc_ + Ak+N) and B = fc_ 

According to [1 2 , in the early stage of aggregation of 
Brownian particles where only monomers (1) and dimers 



FIG. 2: (color online). Time-evolution of the average hydro- 
dynamic radius of the colloidal suspension measured by DLS. 



(2) are present, the temporal evolution of the average 
hydrodynamic radius of the system as measured by dy- 
namic light scattering (DLS) is given by 



I 



h(q) ni (t) 



h{q)n 2 {i)/r ha 
h(q)n 2 {t) 



(4) 



where Ii(q) and hiq) are the wave- vector (g)-dependent 
intensities of radiation scattered by monomers and 
dimers respectively, while r^.i — 85nm and 77^2 — 
1.38 x 85nm are the hydrodynamic radii of monomer and 
dimer respectively [S] . The temporal evolution of can 
be calculated by substituting Eq. (3) together with the 
conservation relation n 2 (£) = (N — ri\{t))/2 into Eq.(4). 
The resulting expression can be expanded and as we are 
interested in the initial kinetic behavior, we can truncate 
to first order in t. This gives 
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(5) 

Upon taking the derivative and rearranging terms we ob- 
tain the standard form 
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with the effective association constant taking account of 
reversibility given by 



\2k- + 4k+N - 



l6klN 2 
k-(h(q)/h(q))(r hil /r h . 2 )]* 



(7) 

This equation is our key result which allows to extract 
the microscopic dissociation rate fc_ from the K e g mea- 
sured in the experiments. One can verify by taking the 
fc_ = limit of no dissociation that Eq. |7| correctly 
recovers the well-known result for the irreversible ag- 
gregation kinetics [T^]. We access K e g experimentally 
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by applying Eq. (6) to data sets of the kind shown in 
Fig. 2. In the analysis, Vh,i, Th,2, Ii{q), and hiq) are all 
known parameters, as well as k + which for attractive col- 
loids is given by the Smoluchowski diffusion-limited rate, 
k + = (8/3)k B T/rj, with -q the water viscosity [ID] . (The 
details of the fitting procedure are given in [8].) The 
effective association constant is plotted in terms of the 
colloidal stability coefficient (TUj, W e s = k + /K c g, as a 
function of T in the left inset of Fig. 3. It is seen that it 
decreases very sharply around T c , after which it reaches 
a diffusion limited-like plateau. 

The knowledge of fc_ also allows us to estimate the 
energy scales for attraction between the nanoparticles, 
as we now explain. The rate at which a dimer disso- 
ciates into two monomers must coincide with the rate 
at which one particle escapes from the attractive poten- 
tial well that binds it to a second particle. The escape 
is a stochastic process promoted by the Brownian mo- 
tion in competition with the attractive interaction |13j . 
In order to keep the treatment analytical we assume an 
effective attractive square-well potential of width 5 and 
depth V m - m . Such an effective potential is schematically 
depicted in Figjl] The approximate Kramers formula for 
the rate of escape from a square- well potential gives [14] 



(8) 



where D = k B T/6ni]rh,i is the diffusion constant (D ~ 
3.4 x \Q- 12 m 2 /s at T ~ 32°C, in our system). To es- 
timate S, note that the attraction between the particles 
is due to the hydrophobic interaction between the PNI- 
PAM microgel grafted layers. As shown in many studies 
in the past, the typical range of the attractive interac- 
tion between two organic hydrophobic surfaces in water 
is about lOnm [T^. As this value is mostly indepen- 
dent of the surface chemistry [15] , we take it here as the 
width S of the effective well. Furthermore, this value is 
much larger than the molecular roughness of the PNI- 
PAM shell, which is of the order of lnm [5] [3] , such that 
the attractive hydrophobic force dominates over repul- 
sive entropic protrusion forces which are much shorter- 
ranged [15j . There is now only one free parameter in our 
model, the depth of the effective attraction well, V m i n , 
i.e., the binding energy of the dimers. One-parameter 
fits to the measured values of W c g reveal that V m ia has 
an almost step-like behavior, as shown in the right inset 
of Fig. 3. It increases sharply from zero in a narrow T 
interval and saturates at V^? n = Ylk B T . 

The observed sharp behavior of V m i n can be related to 
the T-dependent behavior of thcrmorcsponsive nanopar- 
ticles. In the T interval where the turnover in V m i n oc- 
curs, the PNIPAM network grafted on the particle sur- 
face undergoes a transition from a hydrophilic (L) to a 
hydrophobic (B) state with increasing T. The transition 
is associated with a volume change from Vl to V B . In the 
thermodynamic limit, this would be a first-order phase 
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FIG. 3: (color online). Main frame: The experimental In- 
dependence of the binding energy (symbols) together with 
the theoretical model given by Eq. (13) (line). Triangles: 
2.5 x 1(T 4 . Squares: 1.25 x 10 -3 wt%. Circles: 2.5 x 1(T 3 
wt%. Insets: left-the one-parameter fit to experimental mea- 
surements of Weft as a function of temperature. Triangles: 
2.5 x 1(T 4 . Squares: 1.25 x 1(T 3 wt%. Circles: 2.5 x 1(T 3 
wt%. Line: theoretical fit varying Vmin as the unique fitting 
parameter. The experiments were done using the methods of 
Ref. [12]; right-the average binding energy as a function of T. 



transition |16] with the volume as the order parameter, 
but now it is smeared due to the finite size of a nanopar- 
ticle [T7] . On a phenomenological level, this transition is 
captured by a partition function 

Y(T, P) = JdV cxp[- 3 (T, P, V)/k B T] , (9) 

where g(T, P, V) plays the role of a Landau free energy, 
accounting for the many states of the nanoparticle re- 
alizing volume V. (We use the T — P ensemble since 
the externally imposed parameters are the temperature 
and the pressure P |16j.) The first order-like transition 
means that g(T,P,V) has sharp minima at Vl and V B , 
which are degenerate at the transition temperature T c . 
The partition function can be thus approximated as 



Y(T,P) = Y L (T,P)+Y B (T,P), 



(10) 



where Y a (T,P) = cxp[-«?(P, P, V)/k B T}. The 

probability that the particle is in the a state is 
p a =Y a /Y . In a narrow temperature range \T— T c \ <C T c , 
the restricted partition functions Y a can be written as 



Y a {T,P) =exp 



1 



k B T 



[H a (T c )-TS a (T c )}\. (11) 



Since the width of the transition is only few Kelvins, 



Eq. (11) applies in the whole temperature range from 



PL = 1 to p B = 1- This allows us to identify H a {T c ) and 
S a (T c ) as the enthalpy and entropy of the nanoparticle 
in the a state away from the transition region. In terms 
of these quantities, the probability of the hydrophobic 
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state is p B = (1 + exp[(AH - TAS)/k B T])~ 1 , where 
AH = Hb — Hl > and AS* = S B — Sl<0 characterize 
the change in the enthalpy and the entropy across the 
transition. At T c , one has Yj, = Y B which implies AH = 
TAS. Note that the probability distribution 



1 



PB 



l + exp[-{T c -T)AS/k B T\ : 



(12) 



has the same form as for a two level system, p B = [1 + 
exp(e/fcsr)] _1 with a T-dependent level splitting e = 
-(T C -T)A5. 

In a dilute solution, p B gives the proportion of hy- 
drophobic nanoparticles. To relate Eq. (12 1 to the ex- 
periment, we have to now combine the above statistical 
mechanics model with some assumptions about the ini- 
tial aggregation. We base these on the following prop- 
erties: in the hydrophilic state, the particles repel each 
other near contact due to combination of steric and hy- 
dration forces, while two particles that are both in the 
hydrophobic state attract each other via the hydropho- 
bic force. In analogy with the kinetic model, we assume 
that collisions between two hydrophobic particles hap- 
pen with the binding energy V£? n , and take all other 
collisions to be of hard sphere type. Starting from all 
nanoparticles in a monomer state, if after a short time 
there were a binary collisions, p B a of them were bind- 
ing. (We always focus on the early stage of aggregation 
in very dilute systems where binary collisions are dom- 
inant.) In the kinetic model, the number of short time 
collisions is just the number of dimers rii at short times; 
in analogy we take a = n-i . Now requiring that the short 
time total binding energies n2P%(T)V£? n and n2U mm (T) 
of the statistical mechanics and kinetic models are the 
same gives 



VminjT) 
mm 



1 + exp[-(T c - T)AS/k B T] 



(13) 



Eq. ( |13[ ) is compared with the experimental data in Fig. 
(3). The comparison confirms that T c = 32 + 273.16-ftT 
and allows for estimating AS ~ —1420 k B . This can 
be interpreted as the interfacial hydration entropy loss 
when the microgel layer turns from hydrophilic into hy- 
drophobic. The corresponding energy is —T C AS = 3.62 x 
10 Jmol . The hydration entropy loss per monomer 
was measured for free PNIPAM chains in water and is 
~ 10JK _1 mol _1 [18]. Thus, if chains were to behave as 
free chains, we would have ~ 10 3 equivalent monomers 
per particle. However, chains in the microgel layer are 
grafted and crosslinkcd. The typical entropy change for 
deformed grafted chains is ~ 10 3 times smaller than for 
free chains [19 . The estimate order of ~ 10 6 monomers 
in the grafted layer agrees with the value ~ 3 • 10 6 based 
on the known composition of our system [5] . 

To conclude, we have presented a kinetic model that 
can account for spontaneous dimer dissociation which, 



combined with light-scattering experiments, provides a 
robust protocol for quantifying the dissociation rate and 
the reversible association (dimerization) energy of mu- 
tually attractive Brownian molecules and nanoparticles. 
We have demonstrated our method on PNIPAM-coated 
core-shell thermosensitive nanoparticles where the at- 
traction stems from the hydrophobic effect. We found 
that the association energy changes sharply from zero 
to Ylk B T across the transition temperature for the 
swelling of the PNIPAM shell. Using a statistical me- 
chanical analysis, we have related this behavior to the 
two-level nature of the nanoparticles. Our method opens 
up the perspective of systematically quantifying the ef- 
fective hydrophobic attraction [20] . between complex 
biomolecules, such as proteins. This has potential appli- 
cations to other problems where dissociation events are 
important, such as the amyloid formation T\, responsible 
for neuro-degenerative diseases [4]. 

We acknowledge financial support by the Swiss Na- 
tional Science Foundation, grant PBEZP2 - 131153 
(A.Z.), by the Adolphe Merkle Foundation (J.C.),by EP- 
SRC Grant EP/F032773/1 (B.B.), and by the Deutsche 
Forschungsgemeinschaft (M.B.). 



[1] A. Stradner, et al. Nature 432, 492 (2004). 

[2] M. Ballauff and Y. Lu, Progr. Polym. Sci. 36, 767 (2011). 

[3] A.M. Morris, M.A. Watzky, R.G. Finke, Biochim. Bio- 

phys. Acta 1794, 375-397 (2009). 
[4] CM. Dobson, Protein Pept. Lett. 13 (2006) 219-227. 
[5] M.C. Daniel and D. Astruc, Chern. Rev. 104, 293 (2004). 
[6] T. Cellmer, et al. Biophys. Biochem. 125, 350-359 (2007). 
[7] T.P.J. Knowles, et al. Science 326, 1533 (2009). 
[8] See Supplementary Material at ... 

[9] J.J. Crassous et al, J. Chem. Phys. 125, 204906 (2006); 
J.J. Crassous et al, J. Chem. Phys. 128, 204902 (2008); 
J. J. Crassous et al, Colloid Polym. Sci. 286, 805 (2008). 
[10] W. B. Russel, D. Saville and W. R. Schowalter, Col- 
loidal Dispersions, Cambridge University Press, Cam- 
bridge (2001). 

[11] A. Fernandez-Nieves et al, Langmuir 17, 1841 (2001); 
M. Rasmusson, A. Routh, and B. Vincent, Langmuir 20, 
3536 (2004). 

[12] H. Holthoff et al, Langmuir 12, 5541 (1996). 

[13] H. A. Kramers, Physica 7, 284 (1940). 

[14] P.A. Smith, G. Petekidis, S.U. Egelhaaf, and W. C. K. 
Poon, Phys. Rev. E 76, 041402 (2007). 

[15] J.N. Israelachvili, Intermolecular and surface forces, Aca- 
demic Press, London (2004); J. Israelachvili and N. Pash- 
ley, Nature 300, 341 (1982). 

[16] R. Kubo, et al., Statistical Mechanics, North-Holland, 
Amsterdam, (1999), Chap. 1 . 

[17] The transition from Vl to Vb is superimposed on a much 
slower linear volume decrease which is unrelated to the 
hydrophobicity, hence it can be neglected in our analysis. 

[18] S. Kunugi, et al. Polym. J. 34, 383 (2002). 

[19] K. Ohno, et al. Macromolecules 40, 723 (2007). 

[20] D. Chandler, Nature 437, 640 (2005); M. Athawale et 



al, Proc. Natl. Acad. Sci. USA 104, 733 (2007). 



